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ABSTRACT 

Cluster number counts can be used to test dark energy models. We investigate dark 
energy candidates which are coupled to dark matter. We analyze the cluster num- 
ber counts dependence on the amount of dark matter coupled to dark energy. Further 
more, we study how dark energy inhomogeneities affect cluster abundances. It is shown 
that increasing the coupling reduces significantly the cluster number counts, and that 
dark energy inhomogeneities increases cluster abundances. Wiggles in cluster number 
counts are shown to be a specific signature of coupled dark energy models. Future ob- 
servations could possibly detect such oscillations and discriminate among the different 
dark energy models. 
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1 INTRODUCTION 

Observational measu r ements f rom Sup ernovae 
JPerlmutter et alJ Il999t iRiess et alj 1200 it 12004). Cos- 
mic Microwave Background Radia tion dSpe reel et all 
l2005t) and large scale structures iZehavi et£dT " l2002l) 



strongly indicate the existence of a dark energy component 
which corresponds to ~ 70% of our Universe energy 
budget and is responsible for its current acceleration. 
The most popular candidates to dark energy are the 
vacuum energy, als o dubbed the cosmological constant 
llCarroll et al] Il992t) , a nd scalar fields also known as 
cosm on or quintessen ce JWettericblEoOll: iRatra fc Peebles! 
Il988t IWettericb] 12003) . The cosmological constant is spa- 
tially homogeneous and its equation of state is always a 
constant. Scalar fields have a time-var ying equation of 
state and are spatially inhomogeneous iSteinhardt et al. I 
Il999l : Ist,einhardt~l l2003'l. Different dark energy candidates 
have distinctive astrophysical and cosmological imprints. 
The later are mainly dependent o n the time-e v olutio n 
of the equation o f stat e ( see e.g. |Zlaj^y_et_alJ (1999); 
ISchuecker et al.l J2003t): [Jassa^ Ba ela fc PadmanabharJ 
J2005I) : iDoran. Karwan fc Wettericbl 120051)) and the "be" 
haviour of its perturbat i ons (s e e e.g. IWang fc Steinhardtl 
Jl998ft : lFerreira fc Jovcd dl99gt) : lMa et al. I lll999l) ). 

The redshift dependence of cluster number counts is a 
promising tool to discriminate among different dark energy 
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models. Several a uthors (see e.g. iMultamaki et all |2003); 
ISolevi et all i2005tl ) have already use it to investigate both 
non-c ouple d quintessence models: S UGRA ferax fc Martini 
Il999l) . RP (iRatra fc Peebleslll988b|). and non-standar d cos- 
mologies: Car dassian models llFreese fc Lewis! l2002t) and 
DGP models fovali et al] 120021) . Other" groups have also 
used cluster number counts together with other observables 
to show how future galaxy cluster surveys would constrain 
cosmological parameters like the amount of dark energy to- 
day or the equation of state paramet er (jL^TJeJlioul |200 
Wang S. et~a1l2005l: iLima fc Hulliool lHaiman et. allhop 



Horcl lou fc BergeH2005T) . The possible effects of dark energy 



inhomogeneitie s on cluster abunda nces was investigated by 
iNunes. da Silva fc Aghaniml i2004) for minimally coupled 
dark energy models. However, so far no one has ever tested 
dark energy models coupled to dark matter using cluster 
number counts. 

Scalar field candidates to dark energy coupled to dark 
matter are strongly motivated by extra-dimensional parti- 
cle physics models. A general feature of these theories, is 
that the size of the extra-dimensions is intimately related 
to a scalar field. The later is coupled to all, or a selec- 
tion of matter fields iDamour et al.|| l990al). depending on 
the high energy physics mode l llCarroll et atlll992l: ICarroU 
119921 iBertolami fc Motalll999l) . A non- minimal coupling of 
the quintesse nce field to dark matter is therefore worth 
investigating l|A mendolJ |200d : iMainini fc Bonometol 12004 



Amendola fc Tocchini-Valent ini 2002. 



A i Farrar fc Peeble 



2004; TA^riendTl^^^TocchiTi^^ileritoi il l2002bt I Wcttcrich 
1993:lTocchini-Valentini fc Amendola l2002l)" lt is then nat- 
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ural to think that due to this coupling, inhomogeneities 
in the dark matter fluid will th en propagate to the 
scalar field, affecting its evolution ijNunes fc Motal l200fit 
iBarrow fe Motal I2003T) . Clearly such effect will become 
even more important when dark matter perturbations be- 
come non-linear. Hence, it is interesting to investigate 
the possibility of a dark energy component which may 
present inhomogeneities at cluster sc ales, during the non- 
linear regime of matter pert urbations Jm aor & Lahav 2005; 
iMota fc van de Bruckl2004j . In fact . IWettericbl <200 ll I2OO2T1 
and in T^Ae^^t^ITll2001fl speculated that highly non- 
linear matter perturbations might indeed affect a scalar 
field even on galactic scales. They found that, at least in 
principle, the quintessence field (or a scalar field) could be 



ies. In Padmanabhan & Choudhurv (2002); Padmanabhan 


(2002); 


Bagla et al 


(2003) and ICaussel 


(2003J) more exotic 



models, based on tachyon fields, have been discussed and 
they argued that the quintessence equation of state is scale- 
dependent. This is indeed a general feature of non-minimally 
coupled scalar fields, whose propertie s depend on the lo- 
cal d e nsity of the region they " li ve in" ( Khourv & Wcltman 
20041: IMota fe Barro^l l2004al: iKhourv fe Weltmanl l2003t 
Clifton et all 120041: IMota fc Barrowl l2004bt iBraxetal] 



20041) ). 



In this paper we investigate the possibility of using 
measurements from cluster number counts to differentiate 
among dark energy models. We study quintessence candi- 
dates coupled to dark matter and analyse the cluster number 
counts dependence on the amount of dark matter coupled 
to dark energy. Adding to that, we also consider the possi- 
bility of dark energy models which present inhomogeneities 
at cluster scales, during the non-linear regime of structure 
formation. We then compare with the more popular models 
where dark energy is homogeneous at those scales and where 
inhomogeneities only occur at horizon scales. We conclude 
the article assessing the possibility of near future galaxy sur- 
veys to discriminate quintessence models coupled to dark 
matter. 



2 COUPLED QUINTESSENCE AND THE 
SPHERICAL COLLAPSE 

We consider a fiat, homogeneous and isotropic background 
universe with scale factor a(t). Since we are interested in the 
matter dominated epoch, when structure formation starts, 
we assume that the universe is filled with cold dark matter 
and a quintessence field ((f)). The equation that describes 
our background universe scale factor is (we set % = c = 1 
throughout the paper): 



3H 2 = 8ttG (p m + P4>) 



(1) 



where H = a/a is the Hubble rate, p^ — i<^> 2 + V(<f>) an d 
V(4>) is the scalar field potential. We assume the potential to 
be a pure exponential function V(<f>) = Vo exp(an(j>), where 
k 2 = 8nG. This is widely used in the literature. With the 
correct choice of the pa rameter a this potential leads to 
a late time acceleration I AjnendoTa[]2j)00|: i Brookfield et af] 



l2006l: iBarreiro et alJl20od ICopeland et al.ll2004h . Since we 
are investigating non-minimally coupled quintessence fields, 
Pm, includes both the dark matter coupled to dark en- 



ergy (pcDm) as well as the non-coupled dark matter (p um ). 
Throughout all the paper we use Q m o = 0.3, Q^o = 0.7, 
and h = 0.65. 

It is important to note that our theory dif- 
fers significantly in one key a spect from the work of 
iKaplinghat fc Raiaramanl i2006l) . where instabilities in the 
matter fluid can occur. In our models, the dark energy 
sector is described by a light scalar field, with a mass 
whic h is at most of order H. T he models investigated by 
IKaplinghat fc Raiaramanl (120061) are such that the mass of 
the scalar field is much larger than H for most of its history. 
This can have significant implications upon the behaviour 
of the dark matter background and the growth of perturba- 
tions which may lead to instabilities. 

In order to calculate cluster abundances we need the 
evolution of the linear matter density contrast (5). We 
describe the evolution of an overdensity up to the non- 
linear regime u s ing th e spherical collapse model (see e.g. 
iPadmanabhan I lll995l) ). The radius of the overdense re- 
gion r and density contrast 5 are related in this case by 
1 + S — pmc/pm = (a/r) 3 , where p mc and p m are the en- 
ergy densities of pressureless matter in the cluster and in 
the background, respectively. 

The energy density of cold dark matter in the back- 
ground and inside the collapsing region are simply given b y 
the following analytical solutions (see e.g. lAmendolal ifeoOO")) 



Pum 

PcDM 
Pumc 
PcDMc 



PO^umO ( — ■ I ( — 

\a,iJ \ a 



(2) 



po fi DMO (_j (_j e (3) 



(1 + 5i)po£lumO (~J ' 

(1 + 5j)poflcDMO (~ 



(4) 



„B(j> c )-B{M 



(5) 



where again the subscripts "um" and "cDM" mean uncou- 
pled matter and coupled dark matter, respectively. Uncou- 
pled matter corresponds to both baryons and uncoupled 
dark matter. The function B(<f>) represents the coupling be- 
tween dark energy and dark matte r. We use the same cou- 
pling as in the m odel discussed in l)Holden fc Wandsll2000l : 
lAmendolalkoOOD . B(<f>) = -C*k0, where C is a constant. 
Since our scalar field only couples to dark matter, this con- 
stant sets the ratio of the strength of the dark-dark inter- 
action with respect to the gravitational interaction; It is 
then clearly not constrained by local experiments or by G/G 
measurements 1 . However, it is constrained by primordial nu- 
cleosynthesis bounds on the quintessence energy density at 
that epoch. Notice that, if the baryons were coupled to the 
scalar field as well, then we would need to consider several 
constraints on the coupling which would arise from a va- 
riety of experiments and observa tions of fifth force effects 
iEllis et alJll989t IWettericbllT995l) . 

The total energy densities in the background and inside 
the cluster are respectively, p m = p um + p C DM and p mc — 



1 iDamour etH] il99o3) derived a constraint for dark matter in- 
teraction with a dilaton based on the age of the Universe. This 
constraint assumes a field with no potential and a nowadays mat- 
ter dominated universe, which is clearly not our case. 
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Pumc + PcDMc which evolve accordingly to 

Pm = — 3-p m + ^rPcDM<^, (6) 
a a<p 

r dB • 

pmc = — O-pmc H rT PcDMcpc ■ (7) 

r d<p c 

The equations of motion for the evolution of the scalar 
field in the background and inside the overdensity are in this 
case (iNunes fc Motall2006l) : 



Q d; dV dB 

= —S-4> 7T — -jrPcDM , 

a dtp dq> 

r- dV dB T 

= — •i — <pc t-j r~ PcDMc H — — 

r d(p c dcpc (h c 



where describes 



inside the 
iMota fe var 



quint essence loss o f 
e.g. M aor fc Lahavl 



(8) 
(9) 

energ y 
:2005); 



It is 
pled to dark 
lution J^nendolal 



cluster 
i de Brucld 120041) ) 

known that the system quintessence cou 
matter ha s a scaling attrac to r so- 
20001: [Holdcn fe Wandd l200d 



lAmendola fe Tocchini-Valentinill2002al) with 



C 2 + Ca + 3 7 
(C + a) 2 



7<A = 



37 2 



C 2 + Ca + 3 7 ' 



(10) 



where 7 and 7^ are the background and the scalar field equa- 
tion of state respectively. This attractor has a power law ex- 
pansi on aocC given by p = ^ (l + 2) l|Copeland et alJ 
2004). The solution leads to a late time acceleration for 
p > 1, that is, for a < 2C for a matter background. The 
value of C can be extracted from Eq. IjlOjl where now we 
also have to take into account the coupled and the uncou- 
pled dark matter. Therefore, 



C 



-a + ■ 



a + 



12Q C 



2fi c 



(11) 



is a good estimate for the value of the coupling in the tracker 
regime. Here O c dm = f2cDM/(ficDM + 

The constra ints from nucl e osynthesis imply that 
p.d>(T„s ) < 0.1 iWetterichl ll99Et iFerreira fe Jovcei Il998t 



Sarkarj I l996l). which t ranslates into a 2 > 4:/0,^(r ns 
(Hold en fc Wanda 12000 ). In all our models we choose a = 

10, and C will be chosen according t o equat ion 1111 . 

F ollowing iMaor fc Lahavl l|2005^ : lMota fc van de Bruckl 

(2004), we study the two extreme limits for the evolution 
of dark energy inside the overdensity. In the first case we 
assume that dark energy is homogeneous, i.e. the value of 
inside the cluster is the same as in the background, with 



( ) (P^c 

\a r J 



■P<t>a 



(12) 



Hence, in this case, dark energy perturbations are not 
present at small scales and so <f> c — <p. In the second limit, 
dark energy is inhomogeneous and collapses with dark mat- 
ter. Thus = and C 7^ <p. In this case perturbations in 
the scalar field are important at cluster scales. 

In order to compute the cluster number counts we 
also need the evolution for the linear density contrast (5l)- 
Which is given by 

5 L = -2H(5 L - /) + / 

+ "2" [pm^L + (1 + Sw^^S^ pj, + Sp^Sw^] , (13) 



where 8$ = Sp^/p,/,, with 
5pw, = 4>6d>-\ — r7"50, 



•■w<\, = (J - ,r,J j -l-^Ls<p + 8^ 



and 



f = G 



where 



G(4>) 



dB 



> + 



dB 



(14) 
(15) 

(16) 

(17) 



This system of equations closes with the equation of motion 
for the scalar field perturbations 



= -3H( 



dB 



G pmSh + (8l — f)4> 



'-V 



dB 

dd 



+ \-rr) G(l 



d 2 B „ 



(18) 



Integrating these equations we are now able to obtain 
the growth factor D(z) — Sl(z)/8(0) and the linearly ex- 
trapolated density threshold above which structures will end 
up collapsing, i.e. 8 c (z) = 5l(z = Z co l)- Here z co \ is the 
redshift at which the radius, r, of the overdensity is zero, 
and is obtained using the spherical infall model. Both of 
these quantities are needed to compute the number of col- 
lapsed structures followi ng the Press-Schechter formalism 
JPress fc Schechterlll974l) . 

In order to understand the cluster number counts de- 
pendence on the amount of dark matter coupled to dark 
energy and the behaviour of dark energy inhomogeneities 
during the non-linear regime, we investigate four different 
models/cases. We have chosen the models parameters in 
such a way as to have limiting cases. These give us a good 
understanding of the physics behind large scale structure 
formation and coupled quintessence models, being at the 
same time viable cosmological models. We clarify here the 
four cases under study. 

• Model A (Homogeneous Dark Energy with a Large 
Amount of Dark Matter Coupled) : 

All the dark matter is coupled to dark energy, S1 c dm = 0.25. 
Only baryons remain uncoupled fi um = Q,t — 0.05. From 
equation 111 IB one has C = 27.4. In this model we consider 
that dark energy does not cluster in overdense regions. Its 
energy density is the same both in the cluster and in the 
background. Thus is the same as in equation 1121 . 

• Model B (Homogeneous Dark Energy with a Small 
Amount of Dark Matter Coupled) : 

Only a small fraction of the dark matter is coupled, 
f2 C DM — 0.05. The rest is uncoupled matter f2 um + fi(, = 0.25. 
From equation 1111 one has C = 139.9. As in case A, we 
consider that dark energy does not cluster in overdense 
regions. Hence, it is a homogeneous component, with the 
same density all over the Universe. 




Figure 1. Evolution of S c with rcdshift. Model A: ^ 0, 
n cDM = 0.25, f2 um = n b = 0.05. Model B: T 7^ 0, H cDM = 0.05, 
n um + O b = 0.25. Model C: T = 0, n cDM = 0.25, C um = O b = 
0.05. Model D: T = 0, C cDM = 0.05, H um + ft 6 = 0.25. The 
ACDM case, solid-line, is also plotted for reference. 



• Model C (Inhomogeneous Dark Energy with a Large 
Amount of Dark Matter Coupled): 

All the dark matter is coupled to dark energy f2 C DM = 0.25, 
only baryons remain uncoupled f2 um = fi;, = 0.05. From 
equation 11 I > . one has C = 27.4. In this case we consider 
that dark energy clusters in overdense regions. Hence, 
r^, = 0, which means that dark energy collapses along with 
dark matter. 

• Model D (Inhomogeneous Dark Energy with a Small 
Amount of Dark Matter Coupled): 

Only a small fraction of the dark matter is coupled S1 c dm = 
0.05. The rest is uncoupled matter fi um + Q b = 0.25. From 
equation 1111) one has C = 139.9. As in case C, we also 
consider the clustering of dark energy in overdense regions, 
therefore F^ = 0. 

In figure we have plotted S c (z) for several dark mat- 
ter/dark energy couplings and for both homogeneous and 
inhomogeneous dark energy models. It is interesting to note 
the wiggles in 8 C , which are a feature of dark energy mod- 
els coupled to dark matter. These wiggles come from the 
oscillations in the dark energy scalar field around the mini- 
mum of the effective potential. When allowing dark energy 
to clump with dark matter = (inhomogeneous mod- 
els), these oscillations are strongly translated to the matter 
fluctuations and hence appear in 5 C . Notice that, in the ho- 
mo geneous scenario oscill ations are still present (see figure 
5 in lNunes fc Motel i2006l) ): nevertheless they are very sup- 
pressed and could not be appreciated in the plot. 



3 PRESS-SCHECHTER FORMALISM 



Press and Schechter JPress fc Schechterlll974T) . using the 
spherical collapse model, provided a formalism to predict 
the number density of collapsed objects. Several groups 
JCovernato et. alJI 20001: iGross et. al.'lll998t IJenkins et.aH 
l20^nsTringeT^t^d — II2OO5I1 found significant deviations be- 
tween Press-Shcecther predictions and N-body simulations. 
A bett er agreement to the simu lat ions is given, fo r insta nce, 
by the lSheth fc Tormenl dl999l) or lJenkins et. alj feOOlft fits 
for n (m). However, the formali sm bv [p res^g Scliechte'rl 
( 1974) and its extensions by iBond et alJ dl99lT) and 



lLacev fc Cold dl99ot) . even if crude, predicts the evolution 
of the mass function of collapsed objects well enough for the 
purpose of this paper: the study of how dark energy inhomo- 
geneities and dark matter-dark energy coupling influences 
cluster number counts. 

The main assumption in the Press and Schechter for- 
malism is gaussianity of the matter density field. When the 
density fluctuation field 5(x) is smoothed with a top hat 
window of radius R, i.e, when averaged in a sufficiently large 
volume V = ^-R 3 around each point, it follows a Gaussian 
distribution: 



p(S L ,R) = 



2na 



(19) 



where a(R) is the rms of linear fluctuations Sl- Both a(R) 
and 5l axe redshift dependent. The volume fraction of points 
with 5l ^ 5c is 



(20) 



which is assumed to be equal to the mass fraction in bounded 
objects with M > ifp m R 3 . 

We are interested in the comoving number density of 
collapsed objects in a mass range. To obtain this we have 
to take the derivative of /, which gives the mass fraction in 
objects with mass between M and M + dM, and also multi- 
ply by jj, which converts the result into number densities. 
Here p is the comoving matter density. Thus, the prediction 
of the Press-Schechter formalism for the comoving number 
density of collapsed objects is: 



n{M)dM 



2 2- f. ±-dM 
M da dM 

'2 ( 8 r . \ diner 



dluM 



exp 



2a 2 



pdM 

M 2 



(21) 

Note that there is a factor of two introduced to recover 
the mean matter density. This factor can be better under- 
stood when t aking into account the 'cloud-in-cloud' struc- 
ture of halos ijBond et alJ l|l99^ 'l. 

There seems to be some confusion in the literature re- 
garding equation 12 II . We would like to stress that the mat- 
ter density p(z) in this equation is the comoving mean mat- 
ter density at a given redshift. In most cases, it is constant 
and is equal to the present mean matter density, but not 
always. This equivalence is no longer true when one general- 
izes the Press-Schechter formalism to coupled quintessence 
models. In this case one has to bear in mind that p varies 
with redshift directly affecting the prediction of the number 
density of coll apsed objects. 

Following IViana fc Liddiel il999f) we take the variance 
in spheres of radius R to be 

- 7 (ii) 



a(R, z) = as 



R 



D(z), 



where D(z) is the growth factor and 



j(R) = (0.3r + 0.2) 



2.92 + lo gl 



R 



Sh^Mpc 



(22) 



(23) 



where F is the shape parameter of the transfer function . 
Note that, although the formalism bv lViana fc Liddiel (Il999t) 
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could be crude for the present day precision cosmology ex- 
periments, it is good enough for the Press-Schechter formal- 
ism and for the purposes of this paper. We are not seeking 
exact solutions nor precise confrontations with the obser- 
vational data, but to understand the influence of inhomo- 
geneities in dark energy and dark energy-dark matter inter- 
action on cluster number c ounts. 
As in Sug ivamal il995l) we use 



T = Q m hexp 



Q b (l + V2h) 



(24) 



because it takes into account the baryon component. 

The Press-Schechter formalism gives us the comoving 
number density of halos, which we want to compare with 
astronomical data. In order to make this comparison easier 
we convert n(m) to a cluster number counts per redshift and 
square degree with mass Mmin above 2 • 10 14 Mq/i _1 , 




Figure 2. Evolution of the ratio 5 c /ag,D with redshift. Model 
A: T ^ 0, H cDM = 0.25, O um = n b = 0.05. Model B: T ^ 0, 
S^cDM = 0.05, n um + U b = 0.25. Model C: Ts = 0, n cDM = 0.25, 



■ 0.05. Model D: T 



o, n c 



■ 0.05, fi un 



dN 

dz 



dQ, 



ldeg 2 



dV 

dzdQ, 



n(M)dM . 



(25) 



The comoving volume element per unit redshift, 
dV/dz — dQr(z) 2 / H(z) (with r(z) being the comoving dis- 
tance), depends strongly on the cosmological parameters 
and, as we will see, on the coupling between dark matter 
and dark energy. Therefore it plays an important role on 
determining the total amount of cluster number counts. 

In figure [5] we plot S c /asD as a function of redshift 
for several case scenarios. We find that all coupled models 
have a ratio S c /asD below that of the ACDM-model. For 
non coupled models this is the only relevant quantity and 
it would have meant to expect larger halo densities than 
the ACDM model. For coupled quintessence models, how- 
ever, one has also to take into account the redshift evolution 
of the comoving matter density, which plays a very impor- 
tant role as we will see in the next section. In fact, p enters 
both linearly and also through a(R(M, p)) in the equation 
for the comoving number density of collapsed objects (see 
Eq. J2D). 

It is interesting to notice that both S c and D acquire os- 
cillations throw equation I13H . Hence, the typical prominent 
wiggles we saw in S c (see fig. for inhomogeneous coupled 
dark energy models, cannot be seen in figure |2] The reason 
being that oscillations in 5 C are exactly compensated by os- 
cillations in the linear growth factor when calculating the 
ratio Sc/agD. 



Figure 3. Comoving background matter density as a function of 
redshift. There is a decrease of density because of the coupling be- 
tween dark matter and dark energy. Increasing the coupling leads 
to a faster decreasing of the comoving density with redshift . Wig- 
gles are a characteristic signature of coupled quintessence models. 
Notice that in this plot non-coupled dark energy models would 
correspond to a constant line equal to one. 



model 



ACDM (fiducial) 0.9 

A (homogeneous, large amount coupled) 0.843 

B (homogeneous, small amount coupled) 0.892 

C (inhomogeneous, large amount coupled) 0.224 

D (inhomogeneous, small amount coupled) 0.695 



4 CLUSTERS NUMBER COUNTS 
DEPENDENCES 

We choose to normalize all models by fixing the number 
density of halos n(m) at redshift zero. This is the normal- 
ization taken bv lNunes. da Silva fc Aghaniml(l2004l) . At red- 
shift zero all models have the same comoving background 
density p and growth factor D. Therefore the only depen- 
dence on n(m) is through <5 c (0)/<78 (see ea I2H . The nor- 
malization is done by adjusting as in each model such that 
<5 C (0)/(T8 is equal to the fiducial (as = 0.9) ACDM case. 
The table of computed as is presented below. 



4.1 Dependence on the coupling between dark 
matter and dark energy 

The coupling between dark matter and dark energy results 
into several signatures which distinguish these models from 
the minimally coupled ones. The first imprint is associated 
to the comoving density. In non-coupled dark energy mod- 
els, as the universe evolves, the mean matter density of the 
universe (p m ), gets diluted by a -3 due to the expansion. In 
order to account for the expansion effect one constructs the 
comoving matter density p — p m a 3 . For models with no cou- 
pling between dark matter and dark energy, p remains con- 
stant. However, this is not the case for coupled quintessence 
models, as can be seen from equation (6). 
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Figure 4. Comoving volume compared to Einstein-de Sitter vol- 
ume for the four study-models. Since the dark energy clustering 
does not affect the background evolution the difference is due only 
to the coupling. Models A and C with all dark matter coupled 
to dark energy have much more volume than models B and D, 
in which only a small fraction of dark matter is coupled. The 
concordance ACDM model is also plotted for comparison. 



In figure[3]we plot the comoving matter density, in units 
of its present value, as a function of redshift. We can see that 
p decreases with redshift. For models with all dark matter 
coupled to dark energy, p is reduced a factor of 3 at redshift 
1. Since n(m) depends linearly on p, the cluster number 
counts are reduced by the same factor. 

Coupling dark matter to dark energy not only 
changes p but also the expansion history of the uni- 
verse through equation Q (see e.g. [ Amendolal l)200ffl and 
lAmendola fc Tocchini-Valentinil l)2002atl for the evolution 
of background quantities). In figure [I] we plot the value 
of dV/dz for our models A, B, C and D referenced to 
the Einstein-de-Sitter universe. The Concordance ACDM 
model is also plotted for comparison. Note that the volume 
element is a background quantity, therefore the clustering 
of dark energy during the non-linear regime of matter per- 
turbations does not affect it at all. It is clear from the fig- 
ure that different possible expansions of the universe are re- 
flected in the comoving volume element evolution with red- 
shift. Models with more dark matter coupled to dark energy 
have higher values of dV/dz. 

Increasing the value of dV/dz directly translates into 
increased cluster number counts. In fact, this effect is go- 
ing in the reverse direction to the previous discussed one, 
i.e., increasing the volume element actually compensates or 
even overtakes the reduction of the number counts due to 
a decrease in the comoving density p. The combination of 
both effects can be more clearly seen in figure where the 
cluster number counts for square degree are plotted for the 
four models. The ACDM model is also plotted for compar- 
ison. In this figure, one can see that coupled quintessence 
models have less number counts than the fiducial ACDM. 
Actually, increasing the amount of dark matter coupled to 
dark energy leads to a decrease in the number counts ob- 
tained. This is due to the different S c (z) /a(M, z) values and 
the decrease of the comoving matter density, which becomes 
more important than the larger accessible volume. 

A peculiarity of models where dark energy is coupled 
to dark matter are the oscillations present at the cluster 
number counts (see fig. • Wiggles are in fact a common 
characteristic of coupled quintessence models. These appear 




Figure 5. Redshift dependence of the number of clusters of M > 
2.10 14 Mq/i — 1 for square degree. All models are normalized to 
have the same number density of halos today. ACDM case is 
also plotted for reference. Note the wiggles which are a feature of 
coupled dark energy models. 

when the scalar field oscillates about the minimum of its ef- 
fective potential. Due to the coupling to dark matter, these 
oscillations are transfered to the dark matter fluid, and in- 
duce a corresponding osc illation in the p C DM a nd p C DMc 
components (see Fig. 7 of ICopeland et ail teOOtfO . Notice, 
however, that these wiggles are related to the quintessence 
potential form and initial conditions for the sc alar field 
jMota fc van de BrudJl2004t iNunes fc Motall2006i) . For in- 
stance, other coupled quintessence models which would have 
an effective potential without a minimum may not present 
such wiggles. Similarly, a different choice of initial condi- 
tions for the scalar field may lead to the case where the field 
did not have time to reach the minimum of its potential to- 
day. Hence, one would not see the oscillations. Nevertheless, 
fluctuations in the cluster number counts if detected would 
likely indicate the existence of a coupling between dark en- 
ergy and dark matter. 

4.2 Dependence on the dark energy 
inhomogeneities 

From figure it is clear that models with clustering dark 
energy (inhomogeneous models) have more number counts 
than their homogeneous counterparts. This can be under- 
stood looking at S c and the ratio 8 c (z) /agD(z) (see figures 
Qand|2J. When dark energy clusters with matter, it acts as 
a negative pressure slowing the growth of structures. Mod- 
els with a linear growth factor increasing slowly have more 
structure in the past because we normalize all cases such 
that we have the same number density of halos today. In 
fact, the density of collapsed objects is very sensitive to the 
linear growth factor and to the critical density 8 c (z). For 
inhomogeneous dark energy models, it turns out to be sig- 
nificantly lower than the fiducial ACDM model (see figure 
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0. This is also the reason for their low as in the normaliza- 
tion table. 

Wiggles are a common feature for both the homoge- 
neous and inhomogenous cases. There are, however, some 
differences between these cases. While in homogeneous mod- 
els wiggles are basically only present in background quanti- 
ties, i.e, pm (see fig-EU; m the inhomogeneous cases, this is 
not so. Due to the clustering of dark energy, wiggles will also 
quite distinctly appear in clustered related quantities like 8 C 
and the linear growth factor 5l (see figureQl. Independently 
of the clustering properties of dark energy, oscillations in 
cluster number counts will apear, in coupled quintessence 
models, due to oscillations in the background density. These 
oscilations are propagated via the Press-Schechter formal- 
ism. Since for a given radius of a top hat filter, the vol- 
ume fraction of the density field points with Sl > S c corre- 
sponds to the mass fraction of collapsed objects with mass 
M > 4/3irR 3 p. The mass fraction is then converted to num- 
ber density through the background density. 



5 DISCRIMINATING MODELS WITH 
FUTURE SURVEYS 

It is important to estimate whether future surveys mesur- 
ing cluster abundances will be able to discriminate among 
different dark energy models. In order to assess such possi- 
bility we test our dark energy model B: Homogeneous dark 
energy component with a small amount of dark matter cou- 
pled. The aim here is not to perform a detailed analysis but 
to get an idea of the potential detectability of the features 
of c oupled quintessen ce mo dels. 

iBahcall fc Bodd J2003f) have used the abundance of 
massive clusters (m > 8 • 10 14 Mq) in the redshift range 
z — 0.5 — 0.8 to constrain the amplitude of fluctuations 
erg within 10% in the ACDM case. Such uncertainty in erg 
comes from the presence of very large errors in cluster num- 
ber counts, which are big enough for different models to 
survive. Moreover, errors in the mass determination of clus- 
ters also significan tly change the expected number counts 
jLima fc Hull2004) . 

In the near future new surveys are planned to specif- 
ically fi nd clusters in the sky. The South Pole Telescope 
(SPT) JSPT Collaboration: Ruhl J., et al.ll2004l) . which is 
currently under construction, will use the Sunayev-Zeldovich 
effect to find clusters and determine their mas ses. Also the 
recent ly proposed Dark Energy Survey (DES) iAnnis et al I 
(2005)) 2 will observe almost the same region of the sky and 
provide redshifts for those clusters. Both surveys will share 
an area of 4000 square degrees in which two thousand clus- 
ters are expected to be found. Such large numbers will allow 
to better test and discriminate dark energy models. The 
expected errors in redshift for the SPT+DES clusters are 
cr z = 0.02 for clusters with z < 1.3 and a z < 0.1 for clusters 
in the redshift range 1.3 < z < 2. Where a z encompasses the 
68% probability for the redshift being in the z ± a z range. 

In order to explore the potential detectability of wiggles 
in the cluster number counts from these future surveys, we 
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Figure 6. Top: Cluster number counts for model B (dashed line) 
along with a ACDM fit for the model (solid line). The fit is 
done by adjusting its and f2 m . Bottom: Cluster number counts in 
redshift bins for the 4000dcg DES+SPT survey. Horizontal bars 
show shot-noise errors for the bins. Solid line shows the wiggles 
as a difference between model B and the fit in an arbitrary scale. 



fit model B with a flat ACDM model by varying as and 
O m . We simply integrate the cluster number counts for all 
4000 deg DES+SPT survey in redhift bins. We choose to 
use bins of width 0.05 for z < 1.3 and 0.1 within the range 
1.3 < z < 2, to be consistent with the expected observational 
errors. The best fit is obtined by minimizing \ 2 



2 

X = 



E 



iVf 



(26) 



'SN 



2 Dark Energy Survey: |http://cosmology.astro.uiuc. edu/DES/ 
http://www.darkenergysurvey.org 



where Nt is the number of clusters in the ith bin and af N 
is its shot noise error, which is the expected error in any 
counting statistics. The best fit corresponds to Q m = 0.27 
and (J 8 = 0.861. Both the best fit and model B are shown in 
the top panel of figure |S| In the bottom panel we plot the 
difference of the binned cluster counts between the model 
and the fit. This difference is also in redshift bins. In the 
bottom panel we also plot in horizontal bars the shot noise 
error for each bin. The continuous line represents the un- 
binned difference between model B and the fit. For clarity, 
this difference is arbitrarily scaled It is plotted only to see 
the correspondence between bins, wiggles and the smoothing 
due to the binning. 

In this section we are interested in having a broad 
idea about the potential detection of the cluster number 
counts oscillations, which are a special feature of coupled 
quintessence models. To answer this, one could ask how sig- 
nificant is the difference between the model and the fit, given 
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the expected errors for the cluster number counts in each 
redshift bin. The minimum \ 2 f° r this realization is 47.7, 
which gives a probability of less than 5% for the wiggles 
being explained by stochastic fluctuations from the best fit 
ACDM model. We also performed a Kolmogorov-Smirnov 
test, which is less sensitive to the tidal parts of the distribu- 
tion. After smoothing the wiggles signature with a Gaussian 
beam of half-width= 0.05 in redshift to simulate the errors, 
the Kolmogorov-Smirnov test gives a probability ~ 6% for 
ACDM being the underlying model. Hence, both the \ 2 test 
as well as the Kolmogorov-Smirnov test seem to indicate 
that future surveys could possibly detect those oscillations 
in the cluster abundances. 



future surveys could possibly detect such wiggles and may 
be able to discriminate among dark energy models. 
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6 CONCLUSIONS 

In this article we have investigated the possibility of using 
cluster number counts to differentiate dark energy models. 
In particular, we have studied quintessence models coupled 
to dark matter. We have also compared dark energy mod- 
els that can present inhomogeneities at cluster scales, with 
models that are homogeneous at those small scales. The aim 
is to better understand the dependence of the cluster num- 
ber counts on the coupling between dark energy and dark 
matter, and on the dark energy inhomogeneities during the 
non-linear regime of matter perturbations. 

We have shown that there is a significant dependence of 
cluster number counts on dark energy inhomogeneities and 
on the amount of dark matter coupled to dark energy. In- 
creasing the coupling between dark matter to dark energy 
reduces the cluster number counts. This effect is due to the 
decrease of the comoving matter density and the distinctive 
evolution of S c /a in time. Dark energy clustering is shown 
to increase cluster number counts by slowing down the for- 
mation of structure. Hence, depending on the amount of 
coupling between dark energy and dark matter and on the 
clustering properties of dark energy, these effects can com- 
bine together or against each other to strongly increase or 
reduce cluster abundances. 

Oscillations in cluster number counts in redshift, are 
found to be a specific signature of models with dark mat- 
ter coupled to dark energy. In homogeneous dark energy 
models, these oscillations are mainly present in background 
quantities, such as p m , while in the inhomogeneous case the 
oscillations also appear in perturbed quantities, like S c . We 
have shown that such fluctuations are propagated to the 
cluster number counts producing this very peculiar cosmo- 
logical imprint. 

Finally, we investigated the possibility of near future 
observations to discriminate among different quintessence 
models coupled to dark matter. As an example, we have cho- 
sen to test a particular model where dark energy is coupled 
to a small amount of dark matter and where dark energy is 
homogeneous at cluster scales. We fit this model to a flat 
AC DAI- model by varying as and fi m and minimizing y 2 . 
When plotting, in redshift bins, the cluster number counts 
for all 4000 deg of the DES+SPT surveys, wiggles still re- 
main above the shot noise for some bins. In fact the null test 
from the \ 2 gives a probability of less than 5% for these wig- 
gles being a stochastic realization of a AC DM model. Hence, 
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